Lateral Shift Makes a Ground-Plane Cloak Detectable 
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We examine the efTectiveness of the ground-plane invisibility cloak generated from quasiconformal 
mapping of electromagnetic space. This cloak without anisotropy will generally lead to a lateral 
shift of the scattered wave, whose value is comparable to the height of the cloaked object, making 
the object detectable. This can be explained by the fact that the corresponding virtual space is 
thinner and wider than it should be. Ray tracing on a concrete model shows that for a bump with 
a maximum height of 0.2 units to be hidden, the lateral shift of a ray with 45° incidence is around 
0.15 units. 
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The invisibility cloaks formed by the coordinate trans- 
formation method are rapidly developing from narrow- 
band to broadband Jj-,10J . For a given coordinate trans- 
formation from "virtual space" into "physical space," the 
homogeneous Maxwell equations retain their form, with 
only the constitutive parameters changing together with 
the field values. This formal invariance gives rise to a 
powerful technique for designing optical devices, one of 
which is the invisibility cloak. Recently there have been 
reports on the applications of quasiconformal mapping to 

[ achieve broadband invisibility by utilizing aground plane 

'. of the perfect electric conductor (PEC) @, HS-ES] • This 
new strategy is based on the assumption that a suffi- 
ciently minimized anisotropy of transformed media can 
be dropped and consequently isotropic media in place of 
anisotropic media can be used to construct a broadband 
cloak |5i] . Subsequent experiments have implemented this 

I model in both microwave [6| and optical frequencies ^ 
[Toj . Most studies have considered this kind of ground- 
plane cloak as being undetectable. However, the physical 
consequence as well as the validity of this quasiconformal 
mapping technique, i.e., neglecting of the anisotropy and 
replacement of anisotropic materials with isotropic ma- 
terials, have not been thoroughly discussed in previous 

\ studies. 

In this Letter we examine the cloaking effectiveness 
of the quasiconformal- mapped ground-plane cloak which 
omits the minimized anisotropy. We demonstrate that 

I the anisotropy of the transformed cloak medium is of fun- 
damental importance in preserving both the phase and 

I energy propagation of electromagnetic waves. Although 
the isotropic ground-plane cloak is able to largely pre- 
serve the phase propagation, the energy propagation is 
subject to a lateral shift which will render the cloak de- 
tectable. We show in a simplified case that by dropping 
the minimized anisotropy, although it can be arbitrar- 
ily small, a lateral shift on the order of the height of 
the object to be hidden will occur. In the model of the 
ground-plane cloak hiding a bump with the maximum 
height being 0.2 units, we show by ray tracing that the 
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FIG. 1: An example of creating an invisibility cloak by coor- 
dinate transformation. In general the transformation is not a 
conformal mapping and the transformed cloak is anisotropic. 



reflected ray of an incident ray with 45° incidence is later- 
ally shifted on the incident side of the cloak by a distance 
of 0.15 units compared to the ideal case. Our theoretical 
analysis shows that this ground-plane cloak is equivalent 
to a free space thinner and wider than the ideal case, 
causing a lateral shift. This laterally shifting phenom- 
ena is similar to the effect of simply putting a ground 
plane above the object to be cloaked. The conclusion of 
a perceptible lateral shift, being inconsistent with previ- 
ous verifications through ra y tr acing Q , electromagnetic 
numerical simulation [3, I8l-ll0l|. and experiments [g, 0- 
llol |. provides a balanced view on the performance of a 
ground-plane cloak. 

In order to understand the fundamental role of 
anisotropy in transformation-based invisibility cloaks, we 
first recall the required forms of the constitutive param- 
eters of the transformation medium. In the transforma- 
tion medium that is transformed from an empty space 
with permittivity eg (or permeability /io), the trans- 
formed permittivity tensor ise — eoJ-J /\J\ (the perme- 
ability /I has similar expression), where J is the Jacobian 
of the transformation. If we require e to be isotropic, 
the Jacobian must be a constant times an orthogonal 
matrix and the transformation corresponds to a rotation 
with scaling (here we do not consider inversion opera- 
tion). However, for a free standing cloak [1], the trans- 
formation does not correspond to a conformal mapping, 
and therefore orthogonality is not preserved by the trans- 
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FIG. 2: (Color online) Left column: Ray tracing diagrams of 
a ray with incident angle of Qi. Right column: A Gaussian 
beam with incident angle of Qi. (a),(b): A PEC ground plane 
in free space with no cloak. (c),(d): A perfect anisotropic 
cloak squeezing space between z = and z = d into space 
between z = Sd and z — d. (e),(f): Replace the anisotropic 
cloak with isotropic medium. In the illustration, 5 — 0.5 and 
6*^ =45°. 



formation. This has a physical significance and we can 
quickly give an interpretation. Consider Fig. [1] where a 
transformation is applied to a virtual empty space and a 
void is created in the center. In the virtual electromag- 
netic space, we can consider the horizontal lines as rays, 
and correspondingly the vertical lines as phase fronts. Af- 
ter transformation, the mesh is distorted. Take the first 
quadrant as an example. In order to fit the void in, the 
horizontal lines in this region are bent clockwise, while 
the vertical lines at the same location are bent counter- 
clockwise. This means the energy propagation is not in 
the same direction as the phase propagation, which is not 
possible in an isotropic medium. It should be pointed out 
that a strict conformal mapping [2] does not suffer from 
this problem, but is restrictive in the choice of boundary 
conditions. 

For the ground-plane cloak, the anisotropy is also im- 
portant. Let us first consider a simplified CcLSC, 3jS shown 
in Fig. [2] Consider a PEC ground plane in free space 
at the origin, as shown in Fig. HJa), where the trajec- 
tory of a ray with incident angle 9i = 45° is plotted. 
Using the transformation method, it is possible to con- 
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FIG. 3: fc-vectors used in phase matching between the inci- 
dent and cloak media. The solid circle represents the incident 
medium and the solid ellipse represents the cloak medium. 
The dotted circle represents the isotropic medium that re- 
places the cloak medium. 



struct an anisotropic but homogeneous layer to cloak the 
entire plane, as shown in Fig. W^c) . In this transforma- 
tion, the virtual electromagnetic space between z' — 
and z' = d in Fig. [D^a) is squeezed into the physical 
space between z = Sd and z = d in Fig. [2jc) {5 = 0.5 
in the illustration). The effect of such a cloak is to have 
the ground plane appear to be at the original plane at 
z — while in reality it has been moved to the plane 
at z = Sd so that there is room to hide objects between 
planes at z = and z = Sd. The transformation used 
to generate this cloak is a one-dimensional compression: 
X = x' , y = y' , and z = (1 — S)z' + Sd, where the un- 
primed variables are the physical space and the primed 
variables are the space before transformation. The re- 
sulting relative permittivity and permeability tensors are 
I = ]! = xxj^+yyj^ + zz{l — S). Note that we can also 
treat it as an extreme case of the previous ground-plane 
cloak model where a curved bump is changed to be a flat 
one here @. 

We will first discuss how the cloak in Fig. [5Jc) works. 
For the cloak to be perfect, we require the wave vec- 
tors kz — k'z/{l — S) and = k'^, such that the total 
phase after a round-trip propagation in the cloak layer 
is the same as when the cloak is not there, i.e. as in 
Fig. ma). Because of increased k^ after the wave enters 
the cloak layer, the wave vector (as opposed to the ray) 
is refracted as if the wave has entered into an optically 
denser medium with a smaller refracted angle. Second, 
we also require the ray to come out at the same point C 
in both Fig. ^c) and (a) , and therefore the Poynting vec- 
tor, or the ray, should follow the path A-B-C in Fig.[2Jc). 
If we take the view of energy or ray propagation, the light 
seems to be refracted into an optically less dense medium 



3 




K w=4.0 >i 



FIG. 4: (Color online) A ray (blue solid line) incident at 45° to 
the normal onto an isotropic ground-plane cloak. The bump 
has a maximum height of 0.2 units The red dotted line 
represents the trajectory of an undistorted ray in the absence 
of the cloak and the bump. The background represents n^. 
Notice the distorted ray has been shifted on top of the cloak 
toward the incident point by 0.15 units. 



with a larger refracted angle. 

This discrepancy between the refracted angles required 
respectively by the phase and energy propagation is 
solved by the anisotropy of the cloak. In the fc-surface 
diagram (Fig. [S]), the incident medium, free space, gives 
a circle of radius uj/c and the cloak medium, a negative 
uniaxial medium, gives an ellipse with the major axis 
along the —k^ axis. The phase matching condition is 
represented by the line at kx = ki^ and determines the 
transmitted wave vector fcf. Note that kt undergoes a 
refraction as if the medium is optically denser, but since 
the medium is anisotropic, the normal of the fc-surface 
determines the direction of power propagation, which is 
denoted by the arrow, < St >■ Thus, the transmitted 
ray is refracted at an angle 0t > 9i, while maintaining 
the phase in the z direction, just as desired. Therefore, 
the original reflected beam in the region of interest is 
recovered, producing a perfect cloak. 

We can use the exact solution of Maxwell's equations 
for an incident Gaussian beam [ll| to illustrate this point. 
Figure ^b) and (d) show a Gaussian beam incident at 
45° to the PEC ground plane and the same beam inci- 
dent onto the cloak, respectively, corresponding to the 
ray tracing diagrams in Fig. [21[a) and (c). The instanta- 
neous E field is plotted. The parameters of the Gaussian 
beam are so chosen that the phase and energy propaga- 
tion can be clearly seen. It can be seen that the phase 
and energy propagation directions are refracted at dif- 
ferent angles in Fig. [Ifd), resulting in the fields being 
identical to Fig. [^^b) above z — d. 

It is interesting to investigate what will happen if we 
use an isotropic medium instead. The desired anisotropic 
medium has two principle refractive indices til and 
riT @, where in our case = 1 and nx — 1/(1 — ^). 
The refractive index of the candidate isotropic medium 
is chosen to be the geometrical average n = ^JnjjWT Q 



which means an infinitesimal area should be preserved 
in the approximate implementation. Therefore, the ra- 
dius of the circular k surface of this isotropic medium is 
between the shorter and longer axes of the elliptical k 
surface of the anisotropic medium, as the dotted circle 
shows in Fig. [3l With the same incident beam, the phase 
matching condition determines the refracted fct . Because 
of isotropy, the transmitted Poynting power < S't > is 
in the same direction as kt- It can be seen that kt is 
very close to fcj, meaning that the emitting ray will pos- 
sess a phase close to the case of a perfect cloak, while 
< St > deviates from < St > with a large angle A0t, 
meaning that the emitting ray will be laterally shifted by 
a perceptible amount. Figure [2je) and (f) clearly show 
the ray tracing diagram and the Gaussian beam incident 
onto this isotropic medium. It can be seen that although 
the emitting ray as a whole has been shifted by a large 
amount, the phase of this emitting ray does not change 
much. 

We have thus demonstrated the fundamental role of 
anisotropy using a relatively large anisotropy of 2. We 
now examine the effects of small anisotropy. In the case 
of Fig. [D^c), if we let 8 be small, the anisotropy a will be 
small. In this case the anisotropy a is homogeneous and 
to further reduce a at one location will cause increase of 
a at another location. In other words, the anisotropy is 
already minimal. According to Ref . 5] , if a is very small, 
we can drop it and replace the anisotropic material by 
isotropic material with refractive index n — ^nhfiT- For 
a ray incident on the perfect cloak with incident angle di , 
the distance between the incident point and the emitting 
point, i.e. the segment AC in Fig.[2][c) is expected to be 
s = 2(itan0i, since it is equivalent to a virtual space with 
height of d. However, in the isotropic case, this distance 
is s = 2(\-S)d , ""V" < 2d(l-(5)tan6'i. Therefore 

y 1— sin^ dijrL^ 

the reflected ray is laterally shifted toward the incident 
point by s — s > 2(5dtan6'i = 0((5(i), which is proportional 
to the height 5d of the object to be hidden, even if the 
thickness of the cloak d is greatly increased. 

We now consider a curved bump similar to the case in 
Ref. [5] , as shown in Fig. H) The cloak is created from a 
quasiconformal mapping, which omits an anisotropy ratio 
of 1.042 |5|]. The ray tracing exercise in Fig. |4] shows 
that there is indeed a lateral shift of the reflected ray 
comparable to the maximum height of the bump. Note 
that the incident point of the ray on the surface of the 
bump is not on the exact peak of the bump, but a little bit 
in front of the peak. This is also caused by the omission 
of anisotropy. However, if we examine the phase of the 
ray, due to a similar reason of two close wave vectors kt 
and kt as shown in Fig. |3] (in the present case these two 
vectors are much closer due to smaller anisotropy), the 
phase propagation can be preserved to a large extent, 
whose distortion is difficult to perceive. 

Now we discuss the origin of this shift in the quasicon- 
formal mapping-based isotropic ground-plane cloak. In 
Fig. [SI we first show the procedure of such quasiconfor- 
mal mapping. The virtual free space that appears as a 
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FIG. 5: (Color online) Procedure of quasiconformal mapping, 
(a) A rectangle of 1 x m in virtual free space, (b) Com- 
pressed rectangle of 1 x M in the intermediate space, (c) An 
anisotropic cloak created from quasiconformal mapping. 



rectangle with dimensions 1 x m (Fig. [Sja)) needs to be 
compressed vertically to fit the rectangle of dimensions 
1 X M in the intermediate space (Fig. EJb)), where AI 
is the conformal module of the profile of the cloak with 
one curved boundary 5] (Fig. [Sjc)). From the interme- 
diate space to the cloak is a typical conformal mapping 
with slipping boundary conditions. The anisotropy of 
a = m/M in the anisotropic cloak Q is inherited from 
the homogeneous anisotropy of a = m/M in the inter- 
mediate space in Fig. [Ifb) that is similar to the case we 
discussed in Fig. [2][c). The assumption of the quasicon- 
formally mapped isotropic cloak is that the minimized 
anisotropy can be dropped by letting a = I and tak- 
ing the refractive index n = ^JriLTiT- We now study the 
consequence of this operation and consider the inverse 
procedure of Fig. [SJ As shown in Fig. ^c) , every mesh 
point of the anisotropic cloak has two orthogonal direc- 
tions, corresponding to principal directions in the com- 
pressed rectangle (Fig. [SJb)). In the ground-plane cloak 
the anisotropy is because in Fig. [5jc), the metric in 
dv direction and the metric in du direction ^guu 

satisfy the relation ^/gZv = ^\j9uu- When we force the 
two different principle refractive indexes and ut to 



be n = ^JUhnT^ it is equivalent to decreasing the met- 
ric along dv direction and increasing the metric along du 
direction by the same factor of i/a- With the inverse 
conformal mapping, in the intermediate space, the met- 
ric along du direction is increased from 1 to ^/aj while 
the metric along dv direction is decreased from a to -^a, 
meaning now the intermediate space is an isotropic space. 
What happens in the virtual free space is that its height 
becomes thinner and its width becomes wider. It should 
be pointed out that due to the slipping boundary con- 
dition in the quasiconformal mapping, the interface of 
the top of the cloak is not perfectly matched to the free 
space above the cloak. Therefore there is some addi- 
tional distortion caused by this mismatch. However, for 
the current case we neglect this distortion. 

Now we consider the virtual space of the example in 
Fig. m Because of the smaller height, the distance be- 
tween the impinging point and the emitting point of the 
incident and reflected rays in the virtual space (similar 
to AC in Fig. [2]) is 2{h/^/a) tan 45°. In order to fit the 
wider width into the physical space, we need to squeeze 
this segment AC by ^/a to 2{h/a) tan 45°. We also need 
to consider that due to the slipping boundary condition, 
the compression of mesh in the middle is denser that 
other places. So the squeezed segment AC will be a lit- 
tle shorter and the final lateral shift caused by dropping 
anisotropy is larger than 2/itan45° — 2{h/a) ta,nA5° « 
0.12 units. In our ray tracing exercise, the final lateral 
shift is about 0.15 units. Note that this shift at 45° in- 
cidence can be reduced by increasing the width of the 
cloak. However, if the cloak is very wide, the incident 
wave with the largest incident angle can have a very long 
trajectory inside the cloak, which will accumulate a per- 
ceptible lateral shift. 

In conclusion, the invisibility performance of the qua- 
siconformal mapped isotropic ground-plane cloak is ex- 
amined. A lateral shift of scattered waves, being compa- 
rable to the height of the hidden object, will result when 
anisotropy is neglected. This is because the correspond- 
ing virtual space is thinner and wider than the ideal case. 
Ray tracing and theoretical analysis confirm our point on 
a concrete example. 
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